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Relativistic heavy quarks on the lattice 
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We investigate feasibility of relativistic approaches to the heavy quark physics on both anisotropic and isotropic 
lattices. Our peturbative calculation reveals that the anisotropic lattice is not theoretically adovantageous over 
the isotropic one to control mqa errors. We instead propose a new relativistic approach to handle heavy quarks 
on the isotropic lattice. 
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1. INTRODUCTION 

Lattice QCD is expected to provide the oppor- 
tunity of precise evaluation for weak matrix ele- 
ments associated with B mesons from first prin- 
ciples. The main obstacle on the way is the sys- 
tematic errors originating from large m^a errors. 
At present this difficulty is avoided in most cases 
by employing the nonrelativistic approaches |^,^ . 
On the other hand, the use of anisotropic lat- 
tice Q has a fascinating feature to allow us to take 
the continuum limit. However, there exists a the- 
oretical concern whether or not mgUs errors could 
revive perturbatively or nonperturbatively even 
after they are removed classically. Unfortunately, 
our one-loop calculation of the quark self energy 
strongly suggests that radiative corrections allow 
the revival of mqa a errors. 

Having found that the anisotropic lattice is not 
theoretically advantageous over the isotropic one, 
we propose a new relativistic way to deal with the 
heavy quarks on the isotropic lattice. Cutoff ef- 
fects in the heavy quark system are discussed by 
applying the on-shcU improvement program|^,P 
to the finite mqa case. We show that a proper 
adjustment of four parameters in the quark ac- 
tion reduces the remaining cutoff effects to be 
0{(aAQCD)^)- We also demonstrate a determina- 
tion of the four parameters at the tree-level from 
the on-shell quark-quark scattering amplitude. 

In this report we present the salient points in 
our worklql. 



2. ANISOTROPIC LATTICE 

In order to obtain a generic form of the quark 
action on the anisotropic lattice, let us make the 
operator analysis according to the on-shell im- 
provement program I^J^. Under the requirement 
of various symmetries on the lattice, we find the 
following set of operators are allowed up to di- 
mension five: 



dim. 3 : 


0^{x) = 


= q{x)q{x), 


(1) 


dim.4 : 


Oia{x) 


= q{x)^oDoq{x), 


(2) 




Oib{x) 


Y.i^{^)liDiq{x), 


(3) 


dim. 5 : 


OUx) 


q{x)Dlq{x), 


(4) 




OM 


= Y.^<i{x)DU{x), 


(5) 




OUx) 


= iJ2M^)<^oiFQiq{x), 


(6) 




OUx) 


= *E»j9(a;)crjjF»jg(a;), 


(7) 




OUx) 


= J2iQi^)hoDo,7iDi]q{x 


,(8) 



where the subscript denotes the time compo- 
nent, while i,j — 1,2,3 space components. Two 
degree of freedom in the eight coefficients are ab- 
sorbed in the renormalization of the quark mass 
Zm and the wave function Zq. We also find O^a 
and Ose are related to other operators by using 
the eq. of motion and hence they are redundant. 
The remaining four parameters have to be tuned 
to remove 0{at^s) discretization errors. Afterall 
we obtain the following expression for a general 
form of the quark action on the anisotropic lat- 
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tice: 

Sq = atal J2x^(^) hoDo + i^Ei7iA + mo 
-atr/2{Dl + 77^,^?) - attgr/4 [cBvE^J^^JF^ 



-ce{1 + v)J2i<^oiFQi}] q{x) 



13 
(9) 



where the Wilson parameter r is taken arbitrary. 

We can determine the four parameters 77, ce 
and cb at the tree-level from the on-shell quark- 
quark scattering amplitude. The improvement 
condition is that the parameters should be cho- 
sen to remove 0{at) cutoff effects from the scat- 
tering amplitude at the on-shell point. This is 
an extension of the previous work0 done on the 
isotropic lattice to the anisotropic one. Compar- 
ing u(j)')jfj,u{p) {p' ^ p) on the lattice with the 
continuum expression, we obtain v = rj = ce = 
Cb = ^- On the other hand, the ly and 77 pa- 
rameters can be also determined from the quark 
propagator by requiring that it should reproduce 
the form, 



Sq{p) 



1 -*E;.7mPm + "^fl 
-|-(no pole terms) + 0{a^ 



(10) 



with proper choice for Zg, Zm, v and 77. We find 

— 1/2 

Zq = 1 — rmQat/2, Z^a = 1 — rrriQat/'i and 
V = r] — 1 dX the tree-level. It should be noticed 
that the Dirac spinor takes the correct relativistic 
form only when eq. ( [lO| ) is satisfied. 

Expanding one-loop contributions to the quark 
self-energy ttiq) in terms of p and 777,0, we 

— 1/2 

obtain the expressions oi Zq , Z,„, v and 77 up 
to 0{g^at). As an example, the expression of Zm 
is given by 

Z^ = [l+g^CE/{16TT^){~3\og{X^al) + A^°y)} 



[1 + atm (-r/2 + g'Cp / {167r^)Al^j) } , ( 



11) 



where A is fictitious gluon mass to regularize the 
infrared divergence. Figure |l| shows ^ dependence 
of Am\ Here we consider r = 1 and r ~ ^ 
cases with 77 = 1, which satisfy the tree- level on- 
shell improvement condition. We observe approx- 
imately linear dependences on ^ for Am\ which 
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Figure 1 . ^ dependence of Am'' in Zj^ . Errors are 
within symbols. 

tells us that 0{g^a) contributions to Zm are effec- 
tively of order g^mag = g'^mat^. This leads us to 
conclude that mas errors are allowed to revive at 
the one-loop level. This is a reasonable conclusion 
in view of the on-shell improvement. As far as 
we know there is no symmetry on the anisotropic 
lattice to prohibit the higher dimensional opera- 
tors multiplied by (77703)". Unless such symme- 
try is uncovered, the theoretical advantage of the 
anisotropic lattice over the isotropic one would 
never be confirmed. 

3. ISOTROPIC LATTICE 

Let us explain our new relativistic approach to 
control mqa errors for the heavy quarks on the 
isotropic lattice. The basic idea is an applica- 
tion of the on-shell improvement program to the 
finite mga case. Here we assume that the lead- 
ing cutoff effects are /i(r77Qa, g^, log a)(aAQCD)* 
(7 > 0), where f, {mga, g'^, log a){aAQGDy > 
/i+i(mQa,g2,loga)(aAQCD)'+^ and 771q > 
Aqcd ■ fi {niQa, g^ , log a) are smooth and continu- 
ous all over the range of 777Qa and have Taylor ex- 
pansions at toqo = with sufficiently large con- 
vergence radii beyond niQa = 1. We show that 
our formulation removes the cutoff effects up to 
/i (mga, g^, log a)aAQCD with the use of the on- 
shell improvement. The remaining cutoff effects 
/2('T7Qa, 5^, log a)(aAQCD)^ could be removed by 
continuum extrapolation or may be small enough 
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to be neglected. 

In Ref.H, we have derived the on-shell im- 
proved action from the "ofT-shell" improved the- 
ory, which has the axis interchange symmetry. 
Here we alternatively derive it without using the 
axis interchange symmetry from the beginning. 
In this case the discussion is exactly the same as 
for the anisotropic case in the previous section. 
The generic quark action at all order of ma is 
given by 

= O''^ J2x^i^) [70-Do + vY^iliD., + mo 

''art/2Dl - ars/2Y,iDf ~ aigcB /iJ^ijCTijFi] 

-aigcE/2J2iCrotFoi] q{x), (12) 

where we are allowed to choose rt = 1 and the 
four parameters z/, r^, ce and cb are to be ad- 
justed. We remark that the parameter rg is not 
redundant, while it is in Ref.[|. As a consequence 
the "di field rotation" required in Ref.||^ is not 
necessary in our formulation. This difference orig- 
inates from the ways to find the redundant op- 
erators: we employ the classical field equation, 
while Ref.|^ uses the isospectral transformation, 
which is originally proposed in Ref. . Although 
it is remarked in Ref.jsj that this transformation 
should be local and invariant under rotations and 
reflections of the lattice, Ref. |^ employs the field 
transformation which breaks the rotational sym- 
metry. If the transformation is restricted to be 
rotationally symmetric, which is actually propor- 
tional to the classical field equation, the derived 
conclusion is exactly the same as ours. 

The leading cutoff effects /o (mga, g^, log a), 
which reflect discretization errors of the oper- 
ators q{x)q{x), q{x)joDaq{x) and q{x)DQq{x), 
are completely absorbed in Zq and Zm- The 
next leading scaling violations stem from other 
terms in eq.(12). Their magnitude are expected 
to be fi^ {mga, ,\oga){aAQCD) for the Wilson 
quark action and ff^ {mga, g^ ,\oga){aAQCD) 
with ff^ [mqa = 0,(7^, log a) = for the (mass- 
less) 0{a) improved SW quark action. This 
implies that once we fix the pole mass from 
some spectral quantity, the cutoff effects in other 
spectral quantities are at most of order aAqcD, 
not (mga)", for the Wilson and the (mass- 
less) 0{a) improved SW quark actions. If the 



four parameters in eq.(|l2|) are properly adjusted, 
the remaining scaling violations are reduced to 

/2°"^(mQa, g2,loga)(aAQCD)' 

The four parameters z/, rg, ce and cb sX the 
tree-level are determined from the on-shell quark- 
quark scattering amplitude in the same way as 
for the anisotropic case. The improvement con- 
dition that the parameters should be chosen to 
reproduce the continuum scattering amplitude 
gives V = sinh(TOp)/TOp, — cosh(mp)/TOp + 
sinh(mp)(rf — l/mp)/mp, ce = and cb — r^. 
Among these parameters, v and can be de- 
termined by demanding that the tree-level quark 
propagator should reproduce the continuum rel- 
ativistic form. As expected, both methods give 
exactly the same values for v and r^. 

4. CONCLUSIONS 

The anisotropic lattice requires the four param- 
eters to be tuned for the 0(a) improvement even 
in the massless case, and allows the revival of 
mQtts errors through one-loop diagrams. As for 
the isotropic lattice our proposed action, which 
also needs the adjustment of the four parameters, 
can reduce the cutoff effects to be ©((oAqcd)^)- 
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